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The measurement of the 2P^ 2 2 



to 25f / j^ 1 transition in muonic hydrogen by Pohl et al. [l[ and 
subsequent analysis has led to the conclusion that the rms radius of the proton differs from the 
accepted (CODATA) @] value by approximately 4%, corresponding to a 4.9a discrepancy. We 
investigate the finite-size effects — in particular the dependence on the shape of the proton electric 
form-factor — relevant to this transition using bound-state QED with nonperturbative, relativistic 
Dirac wave-functions for a wide range of idealised charge-distributions and a parameterization of 
experimental data in order to comment on the extent to which the perturbation-theory analysis 
which leads to the above conclusion can be confirmed. We find no statistically significant dependence 
of this correction on the shape of the proton form-factor. 



PACS numbers: 36.10.Ee,31.30.jr,03.65.Pm,32.10.Fn 



I. INTRODUCTION 

The measurement and subsequent analysis of the 
2-P^ 2 to ZSf^ 1 transition by Pohl et al. fjj concludes 
that the proton rms charge-radius is approximately 4% 
smaller than previously accepted (as per the 2006 CO- 
DATA 1 value 0). If accurate, this would indicate that 
either QED is an incomplete description of the contri- 
butions to the transition, something has been missed, or 
something has been incorrectly calculated. 

Following the work of Refs. 0-0], we focus on one par- 
ticular contribution to this transition, namely the finite- 
size correction to the 2Pi/ 2 -2Si/2 Lamb shift. We report 
on the form-factor dependence of this correction as well 
as the implications for the analysis leading to the proton 
rms charge-radius. Discussions of the form-factor depen- 
dence of finite-size effects have recently been reignited by 
de Riijula and this work should settle claims made in 
that reference. 



II. NUMERICAL METHOD 

The numerical method used here has recently been 
summarised in Carroll et al. Q and will not be repeated 
here. The highly abridged description is that we use 
an effective Dirac equation for a muon (with a reduced 



mass appropriate to the fi-p system). This is expected 
to provide a precise approximation to the two-particle 
Bethe-Saltpeter equation, yielding accurate muon wave- 
functions for the various potentials studied. 

The eigenvalues for each eigenstate can be calculated 
by inserting the various potentials into the effective Dirac 
equation and integrating itcratively to produce the con- 
verged wave-function. Accuracy is controlled by compar- 
ing the converged eigenvalues with those calculated using 
a virial theorem and errors are conservatively found to be 
±0.5 /icV. 



III. PROTON FINITE-SIZE CORRECTIONS 

As a first non-perturbative approximation, the Lamb 
shift in hydrogen is calculated using the point-Coulomb 
and point vacuum polarization potentials 

V(r) = V a + Vvp(r) 
Za 
r 
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where here m e represents the electron mass which arises 
as this is a consideration of the production of a vir- 
tual electron-positron pair. The momentum integration 
variable has been changed to d(q 2 ), perhaps disguising 
the fact that the lower cut-off of q 2 = 4 corresponds to 
q = 2m e — the energy required to produce the pair. 
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These potentials can be modified to account for the 
finite-size of the proton by convoluting the point po- 
tential with the proton charge-distribution. For exam- 
ple, the modification of the Coulomb potential gives the 
Fourier transform of the Coulomb potential in atoms 



V{q) 



-ZaAit 



G E {f) 



q 



(2) 



where we note that as the energy transfer to the pro- 
ton is essentially negligible, q 2 and the invariant, Q 2 , 
are functionally identical. The coordinate-space poten- 
tial can then be written in terms of a three-dimensional 
Fourier transform of Ge{o[ 2 )'- 



as 



(2tt)= 



, . Za 

V c (r) = > -Za 
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Since the potential of Eq. (H|) involves the pro- 
ton charge-distribution p(r) — itself a function of the 
rms charge-radius — this leads to a radius-dependent 
quantity. The dependence on the choice of charge- 
distribution is investigated here. We enforce that a 
charge-distribution must satisfy 



P(r) d 3 



1. 



(5) 



and we can investigate the effect of using different forms 
for p(r). If we use, say, an exponential form (correspond- 
ing to a dipole form- factor, see Eq. (|2T1) ) 



p(r) = Ae 



-Br 



(6) 



then we can define the rms charge-radius via a ratio of 
the moments of the charge-distribution, as 



(r 2 ) 



r 2 p(r) d 3 r 
p(r) d 3 r 



(7) 



for which the value of A in Eq. ([6]) is arbitrary, and 
which can be rearranged such that we arrive at B = 

y/ 12/ {r 2 } (and A = B 3 /8n in order to correctly nor- 
malise the distribution). The normalized exponential 
charge-distribution is then given by 



Pe(t) = ^e-" r ; 



V = Jl2/(r»). 



(8) 



We can, however, perform the same procedure for al- 
ternative charge-distributions and determine the depen- 
dence on this choice. For a Gaussian charge-distribution 



— corresponding to a Gaussian form-factor — the nor- 
malised form is given by 



3/2 

p G (r)=(|) e- 



rf = 3/2(t*). (9) 



Similarly, for a Yukawa charge-distribution — corre- 
sponding to a monopole form-factor — the normalized 
form is given by 



Mr) = [ \~ ) e~ r, " r /r; v" = \/G/(r 2 p )- (10) 



-in 



To ensure that we are considering realistic distribu- 
tions, we also include in our analysis a charge distribu- 
tion extracted from a fit to experimental data of the Sachs 
electric form factor of the proton [ll| Ge(Q 2 ) given in 
that reference by 



G E (Q 2 ) = 



l + q 6 T + q 10 T 2 + q 14 T 3 



1 + q2T + qiT 2 + q S T 3 + (J12T 4 + <7l6T 5 

(11) 

for which the values of qi are given in Table [H and 
for which r = Q 2 /AMp. This parameterization is con- 
strained at 0(Q 2 ) by 

lim G E (Q 2 ) = 1 - Q 2 ^ + 0(Q 4 ) (12) 



to reproduce 



'{r 2 } = 0.878 fm. The charge distribution 



for this form factor is calculated via a Fourier transform 
of Ge- We will herein refer to this distribution as 'Ge 
fitted'. 

We note several efforts 0, Q to include an additional 
'Darwin-Foldy' (or similarly named) contribution to the 
definition of the charge radius beyond that determined 
in Eq. © and used in Eq. ([121). Wc notc tnat tlie 
Darwin-Foldy term was explicitly calculated by Barker 
and Glover [lj| as part of the Breit potential. Together 
with the other terms in the Breit potential this is al- 
ready included in the recoil correction to the Lamb shift 
as calculated by previous authors (e.g. 0]) and as such 
appears in the complete analysis determining r p (Line 17 
of Table 1 of Ref. [l( supplementary). 

The model charge-distributions used in our analysis 
are plotted for comparison in Fig. [1] and we note the 
striking differences between the shapes below r = 0.8 fm. 
These are plotted again in Fig. [5] where they are weighted 
by 4irr 2 (as used in the normalization) and though the 
differences are reduced, they remain non-trivial. 

We can use these charge-distributions (calculated at 
each of four selected rms charge-radii spanning 0.2 fm 
surrounding the values given in Refs. [l| and (J) to cal- 
culate the finite-size Coulomb potential of Eq. (fj| - 
which is plotted in Fig. [3] — and hence the converged 
Dirac wave-functions in response to this, in order to cal- 
culate the eigen-cnergies A a of the 2S , 1 / 2 and 2P 1 / 2 eigen- 
states for each charge-distribution. We can then calculate 
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FIG. 1: Comparison of exponential, Yukawa, Gaussian, and 
Gs-ntted charge-distributions, each normalized to unity as 
per Eq. ((5|, calculated for ^/(r|) = 0.878 fm. Note the strik- 
ing differences between the shapes below r = 0.8 fm. The 
dipole distribution is nearest to the Ge-fitted distribution, as 
expected. 



FIG. 2: Comparison of exponential, Yukawa, Gaussian, and 
Gs-fitted charge-distributions weighted appropriately as they 
contribute to the Lamb shift, each normalized to unity as per 
Eq. ©, calculated for \/{r|) = 0.878 fm. Here the differences 
are perhaps not as striking, but still noticeably different. 



the deviation of the Lamb shifts (5 = X2S — ^2p) calcu- 
lated in the point-Coulomb and finite- Coulomb cases to 
determine the magnitude of the correction 



finite 



3 finite 



Jpoint ■ 



(13) 



We calculate the proton finite-size correction to the 
Lamb shift using the aforementioned effective Dirac equa- 
tion method for several choices of charge-distribution 



TABLE I: Coefficients of polynomial fit to Sachs electric form 
factor data for the proton taken from as used in Eq. (|12|l . 
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14.5187212 
40.88333 
2.90966 
99.999998 

-1.11542229 
0.00004579 
0.003866171 

10.3580447 



(viz exponential, Gaussian, and Yukawa) at several 
separated values of the proton rms charge-radius (viz 

Jijfl = 0.7 fm, 0.84184 fm, 0.8768 fm, 0.9 fm). With 

this information, we calculate a polynomial fit to the data 
of the form 



Afi finitc = «(r 2 >+^ 2 ) 3/2 , 



(14) 



in order to compare with other published data. A dis- 
cussion of the role of finite proton size in vacuum po- 
larization potential is given in Ref. 0, Table L. The 
relevant parameters of our fits are shown in Tabic ILTl and 
we find no significant dependence on the shape of the 
proton charge distribution. 

Moreover, we are able to make a comparison to the per- 
turbative finite-size correction (due to the finite-Coulomb 
potential) to the Lamb shift as referenced in Ref. pj and 
derived in full in Ref. [1J| as 



AE 



finite — 



_2tt (ZanY 
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(15) 



with the caveat that this expression was derived for an ex- 
ponential charge distribution (corresponding to a dipole 
form-factor) and does not equally apply for other distri- 
butions. The derivation of Eq. (fT5"|) assumes that the 
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Schrodingcr wave-function is appropriate, in that the 
value at the origin 



l<M0)f 



(16) 



appears, and as such this expression requires a relativis- 
tic correction -FrcI- Alternatively, by numerically cal- 
culating the converged Dirac wave-functions we require 
no such perturbative correction, and a comparison with 
that given in Eq. (fl5|) is consistent. We find agreement 
between our Dirac calculation with an exponential dis- 
tribution and various evaluations of Eq. (|15[) to within 
0.05%, as detailed in Table HB 

As the fit to the experimental form factor data is per- 
formed at a single value of r p we cannot determine a 
polynomial dependence on this quantity. We can however 
interpolate the shifts from our three models to compare 
at a single value of r p . The result of such a comparison 
is that at r p = 0.878 fm, the contribution to the Lamb 
shift due to the finite size of the proton is given by 



a ^Exponential 
finite 



AE. 



Yukawa 
finite 



a TTiGaussian 
A ^finitc 



AE 



Ge fitted 
finite 



(0.878 fm) = 


-3.9850 meV, 


(17) 


(0.878 fm) = 


-3.9830 meV, 


(18) 


(0.878 fm) = 


-3.9868 meV, 


(19) 


(0.878 fm) = 


-3.9799 meV, 


(20) 



in keeping with our conclusion that the form-factor shape 
is of negligible influence. 

We note Ref. [l2[ in which a choice of electric form- 
factor parameterization is compared to an idealized 
dipole form-factor 



Gd(Q 2 ) = (1 



Va 2 



(21) 



and for which the ratio of the two tends to unity at low 
Q 2 . The ratio remains close to unity up to approximately 
1 (GeV/c) 2 , re-enforcing that a dipole is a suitable pa- 
rameterization of the electric form-factor for the purposes 
of this analysis. 



IV. IMPLICATIONS FOR THE PROTON 
RADIUS 

With the parameterizations of the finite-size contribu- 
tion to the Lamb shift determined, it is possible to infer 
a proton rms charge radius (r 2 ) 1 / 2 by re-analyzing the 
measured transition in muonic hydrogen of Ref. [l|. If 
we take all other contributions to the transition prima 
facie (we note that because of the unknown magnitude 
of off-shell corrections to the photon-nucleon vertex [l5[ 
such an analysis is physically inappropriate) we can solve 
the cubic equation 
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FIG. 3: Comparison of exponential, Yukawa, Gaussian, 
and GE-fitted finite-size Coulomb potentials calculated for 
y/Jrfj = 0.878 fm. 



(where we note that a' and b' may also account for finite- 
size effects in the 25 hyperfine splitting Q, not included 
here, and thus a' — a, b' — b) in which the measured tran- 
sition energy is -^measured = 206.2949 ± 0.0032 meV; the 
sum of theoretical radius-independent contributions to 
the transition are L r _ indop = 209.9779 ± 0.0049 meV (in- 
cluding the full Lamb shift and corrections found in |l|); 



TABLE II: Coefficients of polynomial fits to the finite-size 
correction to the Lamb shift AE^ite (refer to Eq. JTJJ) in 
muonic hydrogen for several choices of charge-distribution 
(calculated using the finite-Coulomb potential), and selected 
published values. All values in this table include a radiative 
correction of —0.0275 (r 2 ,) as per Ref. We note that the 
finite-size effects of the vacuum polarization alter these val- 
ues further, and a discussion of this matter can be found in 
Ref. 0- 



-'measure 



d = L 1 -inde P + a'(r 2 )+6'(r 2 ) 3 / 2 , 



(22) 



Name 


p(r) 


a [oc (r 2 )] 


b [oc <r 2 ) 3 / 2 ] 


exponential 


Ae~ Br 


-5.2276 


0.0351 


Yukawa 


Ae~ Br /r 


-5.2275 


0.0378 


Gaussian 


Ae~ Br2 


-5.2276 


0.0323 


Ref. [1] 




-5.22495 


0.0347 


Ref. [3] 




-5.22456 


0.0346 


Ref. [13] 




-5.2249 


0.0363 
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TABLE III: Proton rms charge-radius U (r|) calculated us- 
ing various charge-distributions. In these calculations, the 
remaining analysis of Ref. [1| is taken prima facie, including 
the radiative (and other) corrections to the finite-size effect. 
The errors in the radii calculated here are dominated by the 
experimental error in //measured- Also shown are the values ob- 
tained in Ref. 0] and the previously accepted 2006 CODATA 
value [j|. 



Name p(r) Vv*p) [& n ] 

exponential Ae~ Br 0.84174(67) 

Yukawa Ae~ Br /r 0.84194(67) 

Gaussian Ae~ Br2 0.84155(67) 

Ref. [1] 0.84184(67) 

CODATA [2] 0.8768(69) 



and the remaining coefficients are taken from Table [II] 
Of the three solutions to Eq. (|22|) . only one is physically 
meaningful. The physically meaningful value of the pro- 
ton rms charge-radius calculated for each of the choices of 
charge-distribution are given in Table IIIII and compared 
in Fig. |U 

We do not calculate a prediction for the proton rms 
charge radius based on the Ge fitted charge distribu- 
tion as the above analysis requires knowledge of the 
polynomial dependence on this quantity, which is absent 
from this aspect of our investigation. Nonetheless, the 
similarity between predictions of the finite-size effect at 
r p = 0.878 fm as given in Eqs. (fl~71l20"f suggest that no 
significant changes would be found. 

V. CONCLUSIONS 

The dependence of the proton rms charge-radius 
extracted from an analysis of the measured transi- 
tion in muonic hydrogen on the choice of proton 
charge-distribution (and thus form-factor) investigated 
here is shown to be of negligible importance, despite 



the wide range of higher-order moments (given that 
(r 4 ) = 5(r 2 ) 2 /2 for exponential; 5(r 2 ) 2 /3 for Gaussian; 
and 10(r 2 ) 2 /3 for Yukawa charge-distributions) and in- 
vestigation of a realistic charge-distribution based on ex- 
perimental data (of the electric Sachs form factor). 

The analysis of Ref. [l| suggests a discrepancy with the 
2006 CODATA rms charge-radius of 3.99%. Using the 
charge-distributions detailed herein, we have calculated 
a discrepancy with the 2006 CODATA value of 4.02% for 
Gaussian; 4.00% for exponential; and 3.98% for Yukawa 
charge-distributions, indicating no significant variation 
based on this choice alone. 

For the purposes of comparison, it is unclear via the 
references of Ref. [l[ which charge-distribution has been 
used to obtain the value found there. We do conclude 
however that no choice of charge-distribution is likely to 
alter the prediction for the proton rms charge-radius in a 
statistically significant manner, given the wide range of 
shapes investigated here, and the here confirmed insen- 
sitivity of the muonic hydrogen Lamb shift to the shape 
of the proton charge distribution. 

Turning our attention to the recent estimates of Ref. @ 
in which the charge-distribution (via the form-factor) de- 
pendence of the results of Ref. [l[ are challenged — we 
propose that the results found herein quantify this de- 
pendence sufficiently and rule out the shape of the proton 
form-factor as a source of the proton radius discrepancy. 
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FIG. 4: Comparison of calculated proton rms charge-radii (with errors) for various charge-distributions. Also shown are the 
values of Pohl et al. from Ref. [1[ and (inset) the 2006 CODATA value, Ref. [2j]. Error bars are shown for all data here, but 
in many cases they are not visible. It should be stressed that the analysis of this data remains reliant on the as-yet unknown 
contribution of off-shell effects [15(. 



